. . ._ Insight into the Role of Instantons and their Zero Modes 
Q>^ ; from Lattice QCD 

J.W. Negele'^ 

1) I '^Center for Theoretical Physics 

Laboratory for Nuclear Science and Department of Physics 
Massachusetts Institute of Technology, Cambridge MA 02139 



tin 
(N 



^ ' Evidence from lattice QCD calculations is presented showing that instantons and their 

CN . associated zero modes play a major role in the physics of light hadrons and the propagation 

(^ . of light quarks in the QCD vacuum. 
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During the quarter of a century I have known Koichi Yazaki as a colleague, collaborator 
and friend, we have shared an abiding interest in many body physics - understanding 

QhI how strongly interacting systems develop their rich and complex properties from their 
fundamental underlying interactions. Hence, it is a special pleasure and privilege to be 
able to take part in this symposium being held in his honor and to talk about insight my 
collaborators and I have gotten from lattice QCD into the role of instantons and their 
associated zero modes in the physics of hadrons. 

c^ \ One of the great advantages of the path integral formulation of quantum mechanics and 

field theory is the possibility of identifying non-perturbatively the stationary configura- 
tions that dominate the action and thereby identify and understand the essential physics 
of complex systems with many degrees of freedom. Thus, the discovery of instantons in 
1975 [1] gave rise to great excitement and optimism that they were the key to understand- 
ing QCD. Indeed, in contrast to other many body systems in which the quanta exchanged 
between interacting fermions can be subsumed into a potential, it appeared that QCD 
was fundamentally different, with topological excitations of the gluon field dominating 
the physics and being responsible for a host of novel and important effects including the 6 
vacuum, the axial anomaly, fermion zero modes, the mass of the rj', and the chiral conden- 
sate. However, despite nearly a quarter of a century of theoretical effort, it has not been 
possible to proceed analytically beyond the dilute instanton gas approximation [2] . In the 
intervening years, the instanton liquid model [3-5] provided a successful phenomenology 
and qualitative physical understanding, but a quantitative exploration of the role of in- 
stantons in nonperturbative QCD has had to wait until lattice QCD became sufficiently 
sophisticated and sufficient resources could be devoted to the study of instanton physics. 
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Our basic strategy will be to reverse the usual analytical process of calculating station- 
ary configurations and approximately summing the fluctuations around them. Rather, 
we will use Monte Carlo sampling of the path integral for QCD on a lattice to identify 
typical paths contributing to the action and then work backwards to identify the smooth 
classical solutions about which these paths are fluctuating. The physical picture that 
arises corresponds closely to the physical arguments and instanton models of Shuryak 
and Diakonov[3-5] in which the zero modes associated with instantons produce localized 
quark states, and quark propagation proceeds primarily by hopping between these states. 

1.1. Aspects of continuum instanton physics 

To put lattice investigations in context, it is useful to recall relevant aspects of contin- 
uum instanton physics. Working in Euclidean time, we evaluate a path integral of the 
form JT)[A]eJ '^ ^'^^ ', Hence, as in statistical mechanics, the weight of a configuration 
depends not only on its energy, but also on its entropy - the number of ways it can be 
realized. In addition, tunneling solutions arise as periodic classical solutions in an inverted 
potential. Hence, we may expect BPST instantons [1] which connect degenerate minima 
of differing winding number with the self-dual gauge potential A",(x) = ^l^"^2 having 
scale-invariant action 5* = | / d^x F^^F^^ = Sn'^g''^ = Sq to have a significant presence 
in the vacuum due to the high entropy associated with translation, size, and color ori- 
entation. The action and topological charge density are localized around the center as 

±^;.^;. = ^;.^;. = jM$r ■ ^^^ ^^^ ~ C-fr^p'x^'^expi-J^UAp)'^, where 
the prefactor and the running of the coupling constant in the last factor produce a dis- 
tribution of instantons ~ p^ for SU(3). Physically, we expect this distribution to be cut 
off at large p by interactions between instantons and by fluctuations when the amplitude 
of a sufficiently large instanton becomes small relative to quantum fluctuations. It is the 
difficulty in treating these infrared effects that has stymied analytic progress. 

From the axial anomaly, d^J^fi^lfMlbi^ = 2mX]/V'75V' + i^F^^F^,^, the topological 
charge satisfies the index theorem and, for periodic systems, may be expressed in terms 

of fer,mou eigenfuuctions Q = ^ / FF = n. - „„ = mE. -^|S?^. where „, and 
ur denote the number of fermion zero modes. For an isolated instanton, the zero mode 
is ipQ^x) = ^^Ta "t,T3/2 0. In the limit of light quarks, the Greens function for Nj quarks 

reduces to the product of zero modes n/det[-p+ m\il) j{x)il) f{y) -^^ 11/ V^o(^)V^o(l/) ^'^d 
gives rise to the 't Hooft interaction. Thus, light quarks propagate by zero modes which 
in turn arise from instantons. Based on large N arguments, the Veneziano-Witten for- 
mula [6,7] relates the rj' mass to the topological susceptibility in the pure gluon sector 
X = / ^{Q{x)Q{Q)) = ^{ml+ml,-2m\) yielding the expectation that x = (180 Mel/)^ 

The instanton liquid model [3-5] provides an economical phenomenology of instanton 
mediated quark propagation in the QCD vacuum. The integral over all gluon fields that 
one evaluates in lattice QCD using an ensemble of configurations sampling the action is 
replaced by an ensemble of instantons and antiinstantons of size p ~ 1/3 fm and density 
n ~ 1 f m randomly distributed in space and color orientation, where the values of p 
and n are determined from the physical gluon and chiral condensates. One may think of 
the 't Hooft interaction as a vertex that absorbs left-handed particles of each fiavor and 
creates corresponding right-handed particles, and vice versa for antiinstantons. Mesons 



then propagate in the QCD vacuum by the hopping of quark- ant iquark pairs between 
these vertices, and the quahtative features of the channel dependence arises naturally 
with the pion channel being strongly attractive, the scalar channel repulsive, and the 
interaction in the rho channel very weak. The chiral condensate arises naturally in this 
picture by the fact that the zero modes for isolated instantons mix in the instanton liquid 
giving rise to a finite density of states at low virtuality. 

1.2. Lattice QCD 

A QCD observable is evaluated by defining quark and gluon variables on the sites and 
links of a space-time lattice, writing a Euclidean path integral of the generic form [8] 

= Z-1 /r'(f/)e'"'^^*^(^)'^(^)M-^(f/)M-i(f/) and evaluating the final integral over gluon 
link variables U using the Monte Carlo method. The link variable is U = e*"^^'"^^), the 
Wilson gluon action is S{U) = ^ Y^d{^~^ R-6 Tr Un) where f/n denotes the product of link 
variables around a single plaquette, and M{U) denotes the discrete Wilson approximation 
to the inverse propagator M{U) ^ m +(^ + igjjL. 

Vacuum correlation functions for space-like separated hadron currents calculated in 
lattice QCD display the qualitative behavior expected from the 't Hooft interaction and 
agree semi- quantitatively with the instanton liquid model. As emphasized in ref [9], 
correlation functions of the form R{x) = (0|TJ^(a;) J^(0)|0) characterize the spatial and 
channel dependence of the interaction between quarks and antiquarks and thus supple- 
ment hadron bound state properties like phase shifts supplement deuteron properties in 
characterizing the nuclear interaction. The ratio R{x)/Rq{x) of the interacting to free 
correlator has been calculated in quenched QCD for the following meson and baryon cur- 
rents [10], J = u'yi_,d,u'jf,'j5d,u'j5d,ud,eabc[c''C'y^u'']'yi_,'y5d'', and eabc[u''C'y^u^]u''. Resuhs 
are consistent with dispersion analysis of e"'"-e~ and other data in relevant channels, and 
agree in detail with the channel dependence expected from the 't Hooft interaction and 
the instanton liquid model. Note for subsequent reference, that all correlation functions 
may be decomposed into a continuum contribution concentrated near the origin and a 
resonance contribution arising from the lowest bound state or resonance which dominates 
in the region of 1 fm. Quenched calculations aX [3 = 6.2 [11] corroborate the original 
/9 = 5.7 results. 

2. Extraction of the instanton content of lattice gluon configurations 

2.1. Identifying instantons by cooling 

The Feynman path integral for a quantum mechanical problem with degenerate minima 
is dominated by paths that fiuctuate around stationary solutions to the classical Euclidean 
action connecting these minima. In the case of the double well potential, a typical Feyn- 
man path is composed of segments fiuctuating around the left and right minima joined 
by segments crossing the barrier. If one had such a trajectory as an initial condition, 
one could find the nearest stationary solution to the classical action numerically by using 
an iterative local relaxation algorithm. In this method, which has come to be known as 
cooling, one sequentially minimizes the action locally as a function of the coordinate on 
each time slice and iteratively approaches a stationary solution. In the case of the double 



well, the trajectory approaches straight lines in the two minima joined by kinks and anti- 
kinks crossing the barrier and the structure of the trajectory can be characterized by the 
number and positions of the kinks and anti-kinks. 

In QCD, the corresponding classical stationary solutions to the Euclidean action for 
the gauge field connecting degenerate minima of the vacuum are instantons, and the 
analogous cooling technique[12,13] reveals the instantons corresponding to each gauge 
field configuration. 

The results of using 25 cooling steps as a filter to extract the instanton content of a 
typical gluon configuration are shown in Fig. 1 of Ref. [14] using the Wilson action on a 
16^ X 24 lattice at Q/ g"^ = 5.7. As one can see, there is no recognizable structure before 
cooling. Large, short wavelength fiuctuations of the order of the lattice spacing dominate 
both the action and topological charge density. After 25 cooling steps, three instantons 
and two anti-instantons can be identified clearly. The action density peaks are completely 
correlated in position and shape with the topological charge density peaks for instantons 
and with the topological charge density valleys for anti-instantons. Note that both the 
action and topological charge densities are reduced by more than two orders of magnitude 
so that the fiuctuations removed by cooling are several orders of magnitude larger than 
the topological excitations that are retained. 

2.2. The QCD Vacuum 

Several extensive studies of the instanton content of the SU(3) vacuum have been 
performed recently and are reviewed in ref [15]. The results consistently indicate that 
the average size of an instanton when extrapolated to the uncooled vacuum is 0.39 ±.05 
fm, in agreement with the liquid instanton model, and that the topological susceptibility 
is X ~180 MeV, in agreement with the Veneziano-Witten formula. 

2.3. Comparison of results with all gluons and with only instantons 

One dramatic indication of the role of instantons in light hadrons is to compare observables 
calculated using all gluon contributions with those obtained using only the instantons 
remaining after cooling. Note that there are truly dramatic differences in the gluon 
content before and after cooling. Not only has the action density decreased by two 
orders of magnitude, but also the string tension has decreased to 27% of its original value 
and the Coulombic and magnetic hyperfine components of the quark-quark potential are 
essentially zero. Hence, for example, the energies and wave functions of charmed and B 
mesons would be drastically changed. 

As shown in Fig. 1, however, when the coupling constant, or equivalently, the lattice 
spacing, and quark mass are set by the nucleon and pion masses in the usual way, the 
properties of the rho meson are virtually unchanged. The vacuum correlation function in 
the rho (vector) channel and the spatial distribution of the quarks in the rho ground state, 
given by the ground state density-density correlation function[16] (p|^7o?(a;)g7o5'(0)|p), are 
statistically indistinguishable before and after cooling. Also, as shown in Ref. [14], the rho 
mass is unchanged within its 10% statistical error. In addition, the pseudoscalar, nucleon, 
and delta vacuum correlation functions and nucleon and pion density-density correlation 
functions are also qualitatively unchanged after cooling, except for the removal of the 
small Coulomb induced cusp at the origin of the pion. Similarly, the axial charge matrix 
elements specifying the spin content of the nucleon, {PS\{qYi'y5q)\PS) = 2S''Aq, are 
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Figure 1. Comparison of rho observables calculated with all gluon configurations and 
only instantons. The upper left-hand plot shows the vacuum correlator in the rho chan- 
nel calculated with all gluons and the upper right-hand plot shows the analogous result 
with only instantons. The lower plot shows the ground state density-density correlation 
function for the rho with all gluons (solid circles) and with only instantons (open circles). 

quite similar when calculated with all gluons and only instantons [17]. 



3. Quark Zero Modes and their contributions to light hadrons 

3.1. Eigenmodes of the Dirac operator 

In the continuum limit, the Dirac operator for Wilson fermions, 
Dipx = ^x-fiY.^l[{r-J^,)u^,^^px+lM+{r+-f^)ul,_^f,^f,^p^_^], approaches the famihar continuum 

result ^ m + i{;^ + gj^L) ip . 

In the free case, the continuum spectrum is — [m + i|j}|] and the Wilson lattice operator 
approximates this spectrum in the physical regime and pushes the unphysical fermion 
modes to very large (real) masses. In the presence of an instanton of size p at x = 0, it is 
shown in Ref. [18] that the lattice operator produces a mode with a real eigenvalue which 
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Figure 2. Contributions of low Dirac eigenmodes to the vector (upper graph) and pseu- 
doscalar (lower graph) vacuum correlation functions. The upper graph shows the contri- 
butions of 16, 32, 64, 96, and 128 eigenmodes compared with the full correlation function 
for an unquenched configuration and the lower graph compares 128 eigenmodes with the 
full correlation function for a quenched configuration. 



modes goes to zero as the lattice volume goes to infinity. In addition, instanton-anti- 
instanton pairs that interact sufficiently form complex conjugate pairs of eigenvalues that 
move slightly off the real axis. Thus, by observing the Dirac spectrum for a lattice gluon 
configuration containing a collection of instantons and anti-instantons, it is possible to 
identify zero modes directly in the spectrum. 



Ref. [18] shows the lowest 64 complex eigenvalues of the Dirac operator on a 16^ un- 
quenched gluon configuration for 6/(7^ = 5.5 and k = 0.16, both before and after cooling 
(where 100 relaxation steps with a parallel algorithm are comparable to 25 cooling steps). 
The cooled plot has just the structure we expect with a number of isolated instantons 
with modes on the real axis and pairs of interacting instantons slightly off the real axis. 
However, even though the uncooled case also contains fluctuations several orders of mag- 
nitude larger than the instantons, it shows the same structure of isolated instantons and 
interacting pairs, clearly revealing the zero modes in the original, uncooled vacuum. 



3.2. Zero mode expansion 

The Wilson-Dirac operator has the property that D = 'j^D'^'j^, which implies that 
{'^j\l5\'^i) = unless Aj = A* and we may write the spectral representation of the prop- 
agator {xlD^^ly) = J2i —nr:\ — |4llv^ where Aj = A|. A clear indication of the role of zero 
modes in light hadron observables is the degree to which truncation of the expansion to 
the zero mode zone reproduces the result with the complete propagator. 

Fig. 2 shows the result of truncating the vacuum correlation functions for the vector 
and pseudoscalar channels to include only low eigenmodes[18]. On a 16^ lattice, the full 
propagator contains 786,432 modes. The top plot of Fig. 2 shows the result of including 
the lowest 16, 32, 64, 96, and flnally 128 modes. Note that the flrst 64 modes reproduce 
most of the strength in the rho resonance peak, and by the time we include the flrst 128 
modes, all the strength is accounted for. Similarly, the lower plot in Fig. 2 shows that the 
lowest 128 modes also account for the analogous pion contribution to the pseudoscalar 
vacuum correlation function. Similarly, most of the strength of the disconnected graph 
contribution to the rj' correlation function, which should be particularly sensitive to in- 
stantons, is already provided by the lowest 32 eigenmodes, and the fermionic deflnition 
of the topological charge is nearly saturated by the lowest 8 modes [19]. 



3.3. Localization 

Finally, it is interesting to ask whether the lattice zero mode eigenf unctions are localized 
on instantons. This was studied by plotting the quark density distribution for individual 
eigenmodes in the x-z plane for all values of y and t, and comparing with analogous plots 
of the action density. As expected, for a cooled conflguration the eigenmodes correspond 
to linear combinations of localized zero modes at each of the instantons. What is truly 
remarkable, however, is that the eigenfunctions of the uncooled conflgurations also exhibit 
localized peaks at locations at which instantons are identifled by cooling. Thus, in spite of 
the fluctuations several orders of magnitude larger than the instanton flelds themselves, 
the light quarks essentially average out these fluctuations and produce localized peaks 
at the topological excitations. When one analyzes a number of eigenfunctions, one flnds 
that all the instantons remaining after cooling correspond to localized quark fermion 
peaks in some eigenfunctions. However, some fermion peaks are present for the initial 
gluon conflgurations that do not correspond to instantons that survive cooling. These 
presumably correspond to instanton-anti-instanton pairs that were annihilated during 
cooling. 



4. Conclusion 

Altogether, the lattice calculations reported here provide strong evidence that instantons 
play a dominant role in quark propagation in the vacuum and in light hadron structure. 
We have shown that the instanton content of gluon configurations can be extracted by 
cooling, and that the instanton size is consistent with the instanton liquid model and the 
topological susceptibility agrees with the Veneziano-Witten formula. We obtain striking 
agreement between vacuum correlation functions, ground state density-density correlation 
functions, and masses calculated with all gluons and with only instantons. Zero modes 
associated with instantons are clearly evident in the Dirac spectrum, and account for 
the p, TT, and r]i contributions to vacuum correlation functions. Finally, we have observed 
directly quark localization at instantons in uncooled configurations. 
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